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ABSTRACT 


Name of student: Ashish Robert Roll No: 9520602 

Degree for which submitted: M.Tech. Department: MME 

Name of the thesis supervisor: Dr. Dipak Mazumdar, Professor 
Month and year of submission: Jan, 1998 
Text of Abstract begins here 

Visual opacity of molten metals, relatively large size of the industrial metal processing units 
and high operating temperatures often preclude direct experimental observations. As a 
natural consequence of these, mathematical modelling has emerged in recent years as a ke>' 
process analysis tool in the area of metals and materials processing operations. More 
specifically, in the area of steelmaking, transport type models have been frequently applied 
to characterise liquid steel flow in the reactor vessel and to estimate the rates of the 
associated heat and mass transfer processes such as , melting, dissolution, dispersion, 
mixing etc. It is rather well known that various unit processes in steelmaking are inherenth’ 
multi-dimensional ( steady or transient) and multi-phase ( slag-metal-gas) in nature. 
Consequently, in many instances, for effective and realistic process simulation two phase 
three dimensional transport type models are required ( for example, fluid flow in tundish. 
CC mould, degassing units etc.). The present work therefore concerns with the 
development of a steady state, two phase, three dimensional, turbulent flow calculation 
procedure with the primary objective of predicting melt flow and RTD characteristic in 
steelmaking tundish systems. To develop the present two phase model, a combined 
Eulerian-Lagrangian approach has been adopted in which, the trajectories of the individual 
bubles are stochastically determined in space by solving three additional differential 
equations ( essentially Newton's second law of motion). 

Thus as the first step, governing equations of liquid phase and gas phase motions together 
with the equation of tracer dispersion were derived in three dimension in the Cartesian co- 
ordinate system in conjunction with their associated boumiary conditions. To characterise 

i 



liquid phase turbulence, the two equation k-e turbulence model of Launder and Spalding 
was applied. The governing equations flow , turbulence and tracer dispersion were solved 
numerically adapting a control volume based finite difference procedure embodying the 
SIMPLE ( Semi Implicit Method for Pressure Linked Equations) algorithm of Patankar 
and Spalding. Prior to the analysis of melt flow and RTD in tundish, the results from the 
three dimensional model were evaluated against a set of benchmark solutions. This for 
example includes, flow in a cubic cavity, entrance length against Reynolds number in ducts, 
prediction of bubble rise in a stagnant fluid etc. The primary objective of such an exercise 
was to demonstrate the internal consistency of the software developed and to assess the 
accuracy of numerical prediction. Comparisons between bench mark and corresponding 
numerical solutions indicated excellent agreement between the two. Finally, the three 
dimensional turbulent flow model was applied to predict the flow behaviour and RTD 
(Residence Time Distribution) in three different tundish designs. This for example 
included a single strand slab casting tundish, a two strand slab casting tundish and a 
skewed delta shaped six strand tundish, with or without flow control devices ( viz., dams, 
weirs, dams + weirs etc.). For all the tundish designs, numericalh' predicted flow and RTD 
were compared with the available experimental data reported in the literature and 
reasonably good agreement was achieved between the two. In addition to these, the role 
of auxilliary gas injection in the two strand tundish was assessed from a pureh' 
computational stand point with the aid of the present two phase model. 
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CHAPTER 1 


INTRODUCTION TO THE THESIS 


1.1 Introduction 

Physical modelling was traditionally the preferred tool [1] in metallurgy but of late, 
mathematical modelling has become more popular because of availability of several 
software packages, cheaper computational hardware, growing experience with the tackling 
of a broad range of computational problems and better understanding of the underlying 
fundamental concepts. With the development of new materials i.e., ceramics, composites 
etc. and a more challenging market place, the focus of a process metallurgist has undergone 
substantial reorientation and has shifted to cost control and better product quality. These 
warrant better process control and process optimisation . Towards this, mathematical 
modelling, over the last few years has played a key role, as both process control and 
process optimisation require the representation and understanding of the process/operation 
in quantitative form. 

Apart from lower costs, mathematical modelling has other inherent advantages over 
conventional forms of investigation viz., pilot plant scale experimentation and physical 
modelling. They are as follows: 

1. Speed : With the help of a mathematical model an investigator can study various 
configurations with ease, at a remarkable speed . 

2. Information; The data deri’ied from mathematical model are detailed and complete. 
The values of all possible vjiriables in the domain of interest can be calculated and 
this information can be put to much more effective use. Mathematical modelling 
scores over other forms of investigation viz. physical modelling and pilot scale 
experimentation on this aspect. 

3. Simulation flexibility; With the help of a mathematical model, small scale or cold 
models can be studied. Further, various theoretical possibilities such as, constant 



density and infinite reaction rate options can also be evaluated which at times help the 
researcher in developing a better understanding of the subject. 

The general methodology of a mathematical model development is shown schematically 
[2] in Fig.(l.l). This involves identifying the problem, carrying out simple calculations and 
then drawing a plan for detailed experimental and computational work. It needs to be 
stressed here that mathematical modelling together with physical modelling and pilot plant 
investigation needs to be carried out in a complementary fashion. 

Owing to the distinct advantages mentioned in the preceding paragraph mathematical 
model has become very popular among researchers in the field of materials and 
metallurgical engineering. Numerous investigations thus have been carried out and reported 
in the literature by practising and research engineers alike on various aspects of metals and 
materials processing, accompanying benificiation, refining, casting, & solidification etc. 
In the arena of liquid metal processing operations, mathematical modelling has been 
extensively applied to investigate liquid steel flows, alloy motion, alloy dissolution, mixing 
and RTD(Residence time Distribution) etc. This has assumed considerable relevance since 
visual opacity, relatively large size of the processing units aid extremely high operating 
temperature(~1600°) often precludes direct investigations on industrial metal processing 
units. Thus in the present study an attempt has been made to investigate the melt flow and 
Residence Time Distribution characteristics in different designs of steel making tundish 
system, with the aid of a mathematical model, developed by the prent author. Prior to 
outlining the essence of the mathematical model, it however useful to ellucidate the 
characteristic features of the metallurgical tundish system. 

1.2 Characteristics of metallurgical tundish system 

A continuous casting set up essentially consists of three components, the teeming ladle, the 
tundish and the continuous casting mould. The tundish is the intermediate vessel between 
the ladle and the mould which functions as a distributor and a buffer vessel and helps 




Fig. (1.1) : The general methodology of mathematical model development [2] 
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maintain a relatively uniform flow as the teeming ladle is being emptied and replaced by 
a new one. Fig.(1.2) shows a schematic of the continuous casting process. Over the years 
the quality requirements of the final steel product have become very stringent and the 
tundish therefore is considered to offer the last opportunity for the removal and 
modification of impurities and for any adjustment of trim chemical composition by micro- 
alloying. With the emphasis on the novel casting methods like thin slab and strip casting 
in a single step, the role of tundish has become even more critical. 

The performance of a tundish is very closely related to the nature of steel flow as the basic 
tundish design as well as flow modification devices affect liquid steel motion profoundly. 
For example, a well mixed flow in a given tundish would be more suitable for dissolution 
and alloying where as a plug flow is more suited for inclusion separation. Therefore, the 
basic design of the tundish is important towards performing the desired functions 
effectively. To this end, mathematical modelling can give a reasonable quantitative 
understanding of the flows and the associated Residence Time Distribution in the tundish 
and thus, provide useful information on the overall process. This intum can help achieving 
an optimal tundish design. 

When fluid flows through the tundish it is important to establish the time spent in the 
system by individual fluid elements. The "average time" or the "nominal holding time", of 
the fluid in the system is easily calculated from the definition 

^ Volume of fluid in vessel ( 11 ) 

' Volumetric rate of fluid flow 

However, some fluid individual elements may spend a longer, and others a shorter period 
of time in the system. The minimum amount of time a fluid spends inside the vessel is 
called the "breakthrough time". This departure of actual residence times from the mean 
i.e., the distribution of residence times is an important characteristic of the tundish system 
and influences the latter’s performance as a reactor. 


4 




Fig. (1.2): Schematic representation of the continuous casting proces 
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The Residence Time Distribution(RTD) of a fluid flowing through the tundish can be 
determined either experimentally or via mathematical modelling. Experimental 
determination of RTD involves pulse addition of a tracer in the stream entering the vessel 
and simultaneously monitoring the tracer concentration at the exit with respect to time. In 
contrast, RTD can be calculated from the first principles by solving an appropriate species 
conservation equation. The variation of concentration at the exit as a function of time 
deduced either experimentally or theoretically is typically plotted in the form of 
dimensionless concentration vs. dimensionless time. Such a plot is referred to as the " C- 
curve". A typical C curve is shown in Fig.(l .3) . The C curve thus obtained is analysed to 
assess the hydrodynamic performance of a given tundish . The area under the C-curve is 

roo 

unity as all tracer entering the system should eventually leave it / C{t)d-b=: l . From the C- 

Jo 

curve parameters like the minimum breakthrough time ( r„,„, ), the time at which the peak 
concentration is reached ( ), and the average residence time ( ) can be obtained. 
The mean residence time is calculated from the C-curve as: 

_ r 


These can be converted into their equivalent dimensionless forms ( d„i„ = , dp^^k = 

^peak / K j ^mean ^ ^mea„ ^ K) whcrc is thc thcoretical residence time as defined by Eq.(l . 1 ). 
To interpret the results from the C-curve quantitatively the vessel volume can be assumed 
to be made up of three parts namely, a well mixed volume (E„), a dispersed plug volume 
(Vjp) and a dead volume (IQ . The fractional volumes can be calculated as : 






(0 . .0 ,) 

mm peak' 


(1.3) 


(1-4) 


V - 1 - F - F 


(1.5) 




Fig. (1.3) : A typical " C " curve 
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The relative amounts of these volumes determine the efficiency of a given tundish and can 
also be used as appropriate criteria for improving the tundish design or for comparative 
analysis. The preceding formulation assumes that only negligible volume of liquid spends 
a time greater than twice the theoretical residence time in the tundish. 


1.3 Previous work 

The subject of fluid flow and Residence Time Distribution phenomena has received 
considerable attention over the years. Many studies, computational as well as experimental, 
have been carried out in recent years to investigate the influence of design and operating 
parameters on the tundish hydrodynamics and the associated process efficiency. In general, 
two different approaches have been employed by the researcher’s to investigate flow and 
Residence Time Distribution in tundish. These include, 

1 . Physical modelling using either full scale or reduced scale aqueous models 

2. Mathematical modelling 

Numerous physical model studies have been carried out to investigate the hydrodynamics. 
RTD, inclusion separation etc. using aqueous models of steel making tundish system. 
Both isothermal and non-isothermal conditions were studied through water modelling to 
primarily investigate the variation of exit temperature(strand to strand) as a function of 
time. Such aqueous model studies have been carried out either using full scale models or 
appropriately scaled dc wn models. Often the Froude number based modelling criteria ha^■e 
been applied to maintain dynamic similarity between the model and the full scale tundish 
systems. Extensive eicperimental trials were made on flow and RTD embodying various 
kinds of flow modification devices, such as dams, weirs, slotted dams etc. The principal 
objective of the experimental studies was concerned with an optimum tundish design 
together with a set of optimally placed flow control devices. Table 1.1, summarises the 
work reported to date in various aspects of physical modelling of tundisf system. This 



includes the group of researchers together with their area of research in a chronological 
fashion. 


Table 1.1 : Major research efforts in the field of physical modelling of tundish . 


Year 

Researchers 

Area of research 

1981 

F. Kemeny. A. McLean,T. R. 
Meadowcroft and J. D. Young 

Fluid flow studies in a water model tundish of a slab caster wilh/uithoul 
flow control devices such as dams, wiers, and slotted dams. [3] 

1986 

Y. Sahai and R. Ahuja 

Fluid flow and mixing of melt in steel making tundishes with flow 
control devices such as dams , wiers, slotted dams and submerged gas 
injection. [4j 

T. Robertson and A. Perkins 

Physical and mathematical modelling of liquid steel temperature in 
continuous casting tundish. Water modelling experiments were carried 

Ota to develop tundish designs which give optimum temperature over the 
outlet streams of multistrand tundishes.tS] 

1992 

Y. Sahai and 

M. D. Burvval 

Validity of Retuolds and Froude similarity criteria for water modelling 
of melt flow in tundishes. [6] 

Liu Xintian, Zhou Yaohe, Shang 

Baolu, and Jiag V/eiming 

Flow behaviour and filtration of steel melt in continuous casting 
tundish.[7] 

. 1993 

S, Singh and S. C. Koria 

Water model study of the dynamics of flow in a single strand aqueous 
model tundish. [8] 

S. Joo and R. I. L. Outhrie 

Study of inclusion behaviour and heat transfer phenomena in steel 
making tundishes via aqueous modelling. [9] 

199*1 

F. Shen, J. M. Khodadadi, S. J. Pien, 
and X, K. Lan 

Mathematical and physical modelling studies of molten Aluminum flow 
inatundish.[10] 

1995 

Dipak Mazumdar, Guler 

Yanianogiu, Rainani 

Shankamarayanan and Roderick I. L. 
Guthrie 

Similarity considerations in the physical modelling of steel making 
tundish systems illustrating tlie importance of Froude number based 
scaling,[ll] 

1996 

Y. Sahai and Toshihiko Emi 

Sfelt flow characterisation in continuous casting tundishes and 
assessment of mixed flow models of relevance to tundish 
hydrodynamics. [ 1 2] 

1997 

Dipak Mazumdar, Guler 

Yamanoglu, and Roderick I. L. 

Guthrie 

.\ comparative study of hydrodynamic performance of two strand, fiw 
strand, and six strand tundish systems with/without flow control detices 
'ia appropriately scaled down model tundishes.[l 3] 


, Mathematical modelling of melt flow and RTD in steelmaking tundishes were carried out 
through solution of turbulent Navier -Stokes equations in conjunction with an appropriate 
species conservation equation. Mostly bench mark CFD codes such as, PHOENICS. 
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METFLO etc. were applied to this end. 


Table 1.2 : Major research effort in the field of mathematical modelling of tundish. 


Year 

Researchers 

Area of research 

1986 

J. Szekely and N. El-Kaddah 

Mathematical modelling of three-dimensional heat flow , fluid flow and 
nirbulence phenomena in tundishes. [14] 

O. J. Degbusi and J. Szekely 

The modelling of fluid flow, tracer dispersion and inclusion behasiour in 
tundishes.] 15] 

K. Y. M. Lai, M. Salcudean, S, 

Tanaka, and R. I. L. Guthrie 

Mathematical modelling of flows in large tundish systems in 
steelmaking and verification of predicted results svith physical modelling 
results.] 16] 

1987 

K. H. Tacke and J. C. Ludwig 

Study of steel flow' and inclusion separation in continuous casting 
tundishes via mathematical model.] 1 7] 

Y, He and Y, Sahai 

Study of tundish wall inclination on the fluid flow and mixing via water 
and mathematical model. ]1 8] 

1988 

0, J. Degbusi and 

J. Szekely 

Study of fluid flow and tracer dispersion in shallow tundishes via 
mathematical model.] 19] 

1989 

0. J. Degbusi and 

J. Szekely 

Effect of extemaUy imposed magnetic field on tundish performance.]20] 

1991 

0. J, Degbusi and 

J. Szekeh' 

Heat transfer and the role of auxiliary heating in tundishes. ]21] 

1992 

J. L. Yeh, W, H, Hwang, and C. L. 
Chou 

Physical modelling validation of computational fluid dynamics code for 
tundish design. ]22] 

S. Chakraborthy and Y. Sahai 

Mathematical modelling of transport phenomena in continuous casting 
tundishes during ladle transfer operations[23] 

influence of varying ladle stream temperature and]24] 

and of holding time and surface cover in ladles on liquid steel Dow in 
continuous casting tundishes. ]25] 

1993 

S. Joo and R. I. L. Guthrie 

Inclusion behaviour and heat transfer phenomena in steel making 
tundishes via mathematical modelling. [26] 

B. Kaufinann, A. Niedermatr. 

H. Sattler and A. Preuer 

Separation ofnomnetallic particles in tundishes. [27] 

1994 

0. J. Degbusi 

.Application of the Two-fluid model of Turbulence to tundish 
prDblems.[28] 

1993 

C. Damle and Y. Sahai 

The effect of tracer density on melt flow characterisation in continuous 
casting tundishes via mathematical and water model. [29] 


Alternatively Sahai and co-workers developed their custom built model to investigate 
liindish hydrod>Tiamics. In addition to the studies on flow and RTD the role of auxilian" 






































heating and electromagnetic stirring on tundish performance were also studied numerically. 
Some researchers also went on to investigate the floatation characteristic of non-metallic 
inclusion in tundishes. Some recent computational work includes an assessment study of 
the applicability ofk-e model to tundish systems. Mathematical modelling studies carried 
out and reported to date have been summarised in Table 1 .2. This show the various group 
of researchers together with their specific area of research chronologically. 

Alternatively Sahai and co-workers developed their custom built model to investigate 
tundish hydrodynamics. In addition to the studies on flow and RTD the role of auxiliary 
heating and electromagnetic stirring on tundish performance were also studied numerically. 
Some researchers also went on to investigate the floatation characteristic of non-metallic 
inclusion in tundishes. Some recent computational work includes an assessment study of 
the applicability of k-€ model to tundish systems. Mathematical modelling studies carried 
out and reported to date have been summarised in Table 1 .2. This show the various group 
of researchers together with their specific area of research chronologically. 

1.4 Scope of the present work 

The brief introduction to the thesis indicates that mathematical modelling in recent years 
has become very popular to investigate various aspects of liquid steel processing in the 
tundish. Despite numerous modelling exercises taken up over the years in widely varying 
areas, there remains considerable scope in broadening the applicability of the mathematical 
( better called 'Process') models and their performance. A quick survey of the relei'ant 
literature might as well lead one to conclude that mathematical modelling in the area of 
metals processing has been most exhaustive in 2-D or axi-symmetrical systems. Transient 
two phase, 3-D calculations have been few and far between. Since, most metal processing 
operations do not involve symmetry, consequently the need for a three dimensional 
transport model can be readily visualized (for example, analysis of asymmetric gas 
injection, dual plug bubbling, flow in a tundish etc. necessitates a three dimensional 
turbulent flow calculation procedure). 


U 



The purpose of the present work has been primarily to develop steady state, three 
dimensional, two phase, turbulent flow calculation procedure embodying the two equation 
A-'e-model of turbulence in the Cartesian co-ordinate system. The three dimensional code 
developed was extensively tested against available benchmark solutions to assess its 
internal consistency and reliability. Subsequent to these, the 3-D model was applied to 
predict flow phenomena and Residence Time Distribution characteristics in three different 
tundish systems with/without flow control devices. For the first time numerical 
computations were also carried out as exactly as possible to investigate the role of 
auxilliary gas injection in a tundish. Results thus obtained were also compared with those 
published already in the literature. 

1.5 Layout of the thesis 

The main body of the thesis consists of five chapters. Chapter-2 deals with the 
development of the three dimensional, steady state, two phase turbulence flow model. This 
include the governing equations, salient features of the formulation ,boundary conditions 
and the procedure for incorporating inclined walls/flow modification devices in the 
numerical solution scheme. In chapter-3 a rigorous assessment of the computational 
procedure has been carried out by performing some benchmark tests. Subsequent to this 
the model has been applied to predict fluid flow behaviour and Residence Time 
Distribution for different tundish designs with a gas bubbler and various flow modification 
devices. Chapter-4 summarises tlie general conclusions derived from the present work and 
the thesis concludes with the recommendation for future work in Chapter-5. 



CHAPTER 2 


MATHEMATICAL MODELLING 


2.1 Introduction 

The degree of complexity required to mathematically model different metallurgical 
processes varies from one to another. While, operations like axisymmetric gas injection, 
axisymmetric pouring and drainage, can be effectively modelled using a two dimensional 
turbulent flow model, others, like melt flow in tundishes or in a R-H degasser require a 
three dimensional turbulent flow model for effective simulation. Needless to mention that 
transport processes in practically all metallurgical operations can be successfully modelled 
using transient, two phase, three dimensional models. To this end, several software 
packages are currently commercially available. Some of the popularly used software 
packages have been listed in Table-2.1. 


Table-2.1: Popularly used software packages to study transport phenomena [2]. 


Name 

Application 

Technique 

PHOENICS 

fluid dynamics. Heat Transfer, 

Chemical Reaction 

finite volume 

TEACH - 

T 

fluid dynamics 

finite volume 

FLOW3-D 

fluid dynamics 

finite difference 

FLUENT 

fluid dynamics. Heat Transfer 

finite difference 

NECTON 

fluid dynamics. Heat Transfer 

finite element 

NISA 

fluid dynamics. Heat Transfer 

finite dfference 








However, the key problem lies in the definition of the boundary conditions and the abilit}- 
to interpret the results in an intelligent manner. Hence, many researchers prefer custom built 
models from the first principles with consideration of elementary control volumes. Indeed, 
this later philosophy has been applied in the present investigation to develop a steady, Pvo 
phase, three dimensional turbulent flow calculation procedure. 

Thus, the present chapter deals with the development of a steady state, three dimensional, two 
phase, turbulent flow model in Cartesian coordinate system. The computer code developed has 
been subsequently used to study some benchmark conditions together with the fluid flow 
behaviour in different shapes of trmdish with and without flow control devices . 

2.2 Assumptions in modelling 

The following assumptions and idealisations were made in the present work to computationalh’ 
investigate fluid flow and RTD behaviour in steel making tundish systems. 

i. The flow is assumed to be turbulent and macroscopically steady i.e., the effect of 
transience during filling and emptying of tundish was neglected. 

ii. The melt surface is assumed to be flat and phenomena like wave formation, presence 
of slag layer at the free surface and vortexing during drainage etc. are neglected. 

iii. Entrainment of air or gas through the metal inlet is neglected (the ladle stream is 
essentially shrouded ). 

iv. The flow is essentially Newtonian and incompressible and 

V. The system is isothermal and thus the effect of temperature differential on the induced 
flow is neglected. 

In addition to the above, to model the influence of gas injection in tundishes, the following 
assumptions were made; 

vi. Discrete monosize bubbles are assumed to form at the lance/tuyere tip. Furthermore, 
the size of the bubbles forming at the orifice is assumed to be known. 



vii. Bubble-bubble interactions were assumed to be negligible and hence, the drag 
coefficient correlation applicable to a single bubble fluid system is applied for the 
appropriate range of Eotvos and Weber numbers. 

viii. Bubbles within the upwelling gas-liquid plume contribute to the generation of 
turbulence and the turbulence phenomena affects the bubble motion. 

2.3 Description of the steady state, two phase, three dimensional turbulent flow model 

The essence of the combined Eularian-Langrangian two phase, three dimensional turbulent flo^^■ 
model is presented below. 


2.3.1 The liquid phase equations of motion 


In terms of Cartesian co-ordinate system {x,y,z), the governing flow equations ( for the liquid 
phase ) under steady state conditions can be represented as: 

Equation of continuity : 


—{a.^ pu) . pv) - —(a, pw) = 0 
dx dy dz 


( 11 ) 


Equation of motion in x direction : 


-|-(a; pww) * A(a puv) . pinv) = - 
dx dy dz dx 


d f du-. d , du-. d , du^ 

dx dx dv dv dz dz 


(2.2a) 


in which, 


d , d , dv-. d / SWx ^ 


(2.2b) 


Equation of motion in the y direction : 


A(a pjA>) - pw) ^ ^(a, pvw) - - a, 
dx dy dz 


dp 

dy 


d f dv^ 


d f dv-. 


d / dv^ n 

(oc J u — ) + o 

dz^ ^^‘dz^ 


(2.3a) 



in which. 


S = -i.(a n . 


d , dv. 


a f aw,, 


(2.3b) 


Equation of motion in the z direction : 


dx 


(a, pmv) * -!-(« pvw) + -^(a pww) = - ^ }i -^) 

dy dz ‘ ‘ dz dx dx 


d , aw, 
— (a, u — ) 

ay '^"ay 


dz^ ^^^dz^ " 

(2.4a) 


in which, 



d f du-. 


a f av, 

^ dz 


a , aw, „ 


(2.4b) 


Eqs (2.2a) through (2.4b) are the well known Navier-Stokes equations. In these u, v, wand p 
(gauge with reference to local hydrostatic pressure) are respectively the time averaged velocity 
components and pressure, a, is the volume fraction of liquid and F's are the various forces 
acting on the fluid which may include buoyancy, drag etc. Finally, , appearing in Eqs. (2.2a) 
through (2.4b) is the effective viscosity and is derived from an appropriate turbulence model. 


The values of a,, the volume fraction of liquid in a given control volume, in Eqs(2.2a) through 
(2.4b) can be obtained from the solution of bubble trajectory equations (see later), 


2.3,2 The turbulence model 


In the present work, the popular two equation A-e Model of Launder and Spalding[31] has been 
used for the estimation of turbulence parameters. The model incorporates an equation for the 
conservation of turbulence kinetic energy, A { = l/2(u'^ + v'^ - w'^) } and its dissipation rate, 
e{ = -dt’dt }. The A and 6 transport equations in terms of Cartesian co-ordinate system can be 
expressed as following: 


Conservation equation for the turbulence kinetic energy, A: 



^ (p/ “z wfc) - ^(Pz ttz Uj wk) = -^(a, i^-^) 

cy dz dx a. dx 


dx 


d , P. dL 
ay Ofc ^ 


3 , Pe 
oz 0^ dz 


(2.5a) 


in which. 


- P^) 


(2.5b) 


Conservation equation for the dissipation rate of turbulence kinetic energy, e : 


^(Pz «z“«)^ J-(P/ «z^e) 

ax ay 


+ — (p, a,vi'e) = ~(a, 
dz ' ^ dx ^ 


Pc de^ d /■ Pc de^ d ^ 

_ ) . (a _ ) , (a 

a^dx dy o^dy dz 


Pc de 


o„ dz 

5 


)^‘^e 


(2.6a) 


in which. 


a 


C^eG^ - C,pe- 


•] *iPb^l k) “z 


(2.6b) 


Gk appearing in Eqs.(2,5) and (2.6) is the volumetric rate of turbulence generation and for a three 
dimensional Cartesian co-ordinate system, can be expressed as: 




2Pr[(— ) 
uX 


^ dz dy 


dvG ^ rdv ^ dwG ^ 
dx *^dz * dy^ "^dz 


dwG 


(2.7) 


The effective viscosity , p, , appearing in Eqs.(2.2 ) to (2.7) is defined diS,\ = Pt+ pi , in 
which, 



(2.8) 


The five empirical constants of the it-emodel viz,, C/, C,, C^, a^. and <7^ were assigned to their 
standard values [31] . These are summarised in Table-2.2 . As seen fi'om the above, flow and 
turbulence model equations are mutually coupled. 



Table~2.2:VaIues of constant used in a k-e model. 


C; 





1.43 

1.92 

1.00 

1.30 

0.09 


2.3.3 Equation of motion of the gas phase 

The effect of gas injection in the tundish and hence the motion of the gas bubbles in the bulk 
phase liquid has been carried out via a Langrangian approach, in which the trajectories of the 
individual bubbles are stochastically determined in space by solving three differential equations 
in time. The appropriate form of Newton's second law of motion that describes the bubble 
motion, can be represented along the three relevant coordinates, x ( longitudinal), y ( vertical i 
and z ( transverse) according to: 

Longitudinal direction, x; 


(1 . o. 5A)^ . . (Pi, - V ) . 


p; dt 




dy dx dz 


(2.9) 


Vertical direction, y .• 


(1 . 0.5^)^ . . I Pi, - V) . 1.54ilP^ . - W ^] . (1 - 


Pg 


p/" (2.10) 


and transverse direction, z: 


(1 . . af . pp'f ] (2. 1 1) 




The three corresponding kinematic relationships, which define the trajectories of bubbles 











within the two dimensional flow domain are; 


dx 

~dt 



(2.13) 


il-v 

dt ■ 


( 2 . 12 ) 


and. 


^-V 
dt ' 


(2.14) 


The instantaneous drag coefficient, C^, in the Eqs.(2.9) to (2.11) was evaluated using an 
experimentally determined correlation[32] for oblate spheroid / spherical capped bubbles: 


— ; (500 < Re < 5000 ) 

[— 1 0.235] 


(2.15) 


Furthermore, Re in Eqs.(2.9) through (2. 1 1) is the Reynolds number, and is estimated on the 
basis of absolute magnitude of the relative velocity vector, i.e. 


Re = p, d^ [ 




(2.16) 


Similarly in Eqs.(2.9) through (2. 1 1) , U, V and W are the instantaneous components of fluid 
motion, and were estimated by summing up the relevant mean velocity components ( u, v and 
w) with the corresponding fluctuating velocity components assuming turbulent fluctuations are 
essentially isotropic and obey a Gaussian probability distribution function. For example. 



(2.17) 


3 

where i|; is a normally distributed random variable and varies between -1 and 1 . Thus during the 
rise of the bubble through the liquid, the velocity fluctuations were estimated stochastically and 
were assumed to remain constant for a time period, equivalent to the eddy life time, r . The life 
time of an eddy, during which a characteristic velocity fluctuation was assumed to persist, is 
estimated from the turbulence theory according to[32] : 

r = 0.195 A (2.18) 


2.3.4 Estimation of gas volume fraction 

Figure(2.1) illustrates a typical section of the grid network and a bubble trajectory intercepting 
the Eularian grid system. Once the bubble population in the flow domain is predicted from 
solution of the bubble sub-model equations, gas voidages can be readily estimated from the 
principle of volume continuity, according to ; 

a , volume of gas ^ Q,(N * 1)^^ 

^ volume of control volume 2A \ ■ ) 

Further more, liquid volume fraction, a, , can be easily determined from the above since, 

+ (T, = 1.0 . 

2.3.5 Estimation of the interphase drag and the turbulence production due to bubbles 

The volumetric flow rate of turbulence generation via the bubbles, , embodied in the k-e 
model can be represented as : 


Pb = ^b-^ 




4d 


(2.20) 


where, is the resultant bubble velocity and Uj is the resultant fluid velocity. C* inEq.2.20 




@ Scalar nodes O bubbles 


Fig. (2.1) : Schematic of a section of Eulerian grid network intercepted by 
a bubble trajectory ( 2D - representation ). 


is a constant and takes a value between 0 and 1 . The value of Q depends on the physical 
characteristic of the two phase plume ( e.g., bubble population, bubble size etc.). The constant 
Q is normally determined by trial and error. 

Similarly, the momentum exchange terms Fy, and F^ can be evaluated assuming that drag 
forces experienced by bubbles act with equal magnitude but opposite directions in the liquid 
phase. The expressions for F ^ , F^ , and F, can be represented in the compact tensorial 
notation as ; 



q 






( 2 . 21 ) 


2.3.6 Equation of tracer dispersion 

As pointed out previously the residence time distribution (RTD) characteristic of a tundish can 
be studied mathematically by solving an appropriate species conservation equation. Thus, 
following a pulse addition the appropriate conservation equation of the added tracer can be 
described in Cartesian co-ordinate system in terms of the following transient convection- 
diffusion equation; 


dt 


(Q . {-{uQ . |-(vC,) . |-(>vQ 


dx 


dy 


dz 


d 

^(D __I) 
dx ‘ dx ■ 


d ^1 
dy^ ‘‘dy’ 


A(D.^) 

dz ‘ dz 


( 2 . 22 ) 


in which, is the effective diffusivity defined as; = Di + Dr , where D, is the molecular 
diffusivity and Dr is the eddy diflhsivity. The latter, can be directly equated to the turbulence 
kinematic viscosity assuming turbulent Schmidt number to be approximately unity 
i.e., Hr! {pDj) = 1 (valid for a wi.de range of engineering problems). 

2.4 The boundary conditions 

Fig.(2.2) shows the schematic of a tundish for mathematical simulation in the present study. 
The set of boundary conditions applied on the governing equations with relevance to the given 




tundish can be represented as: 


(i) At the free surface ( in the north direction ) zero shear is assumed to be transmitted for the 
velocity components. Zero gradient of ife and ewere specified at the free surface. In conjunction 
with these, vertical component of flow is set to zero. 

Mathematically, 

M.y = H , 0 s j: < L/2 and 0 < - < W/2, 

v = 0, duldy=dwldy = Q , dkldy — deldy^^ 

(ii) No slip condition for the velocity components was specified at the solid walls. The values 
of k and £ were set at zero at the wall. 

At the bottom wall, y = 0 , 0 ^ x < L'2 and 0 < r ^ W/2, 
u-v — w = Q,k-e = 0. 

At the side wall, x = 0. Os ysH and 0 < r s W/2, 
w = v = w = 0 , A:=€ = 0. 

At the frontal wall, r = 0 , 0 < y <R and 0 s x s L/2, 
iy = v = w = 0 , k = e = 0. 

(iii) Gradients of all the velocity components at the symmetry plane are taken as zero and 
together with zero gradient of k and e 

At the central longitudinal symmetrical plane, r = W/2 , 0 s y < H and 0 < x s L/2, 
duldx = dvldy = ^yvIdz-Q , dkicy- d ddy — Q 

At the central transverse symmetrical plane, x = L/2 , 0 < y s H and 0 s r < W/2, 
duldx = dvldy= , dkjcy-deldy=^^ 

|n addition to these, computations were carried out by force fitting the following conditions in 
the nuntorical solution scheme. 



(iv) Inside the pouring stream, 

w = w = 0 , V = 

(v) At the exit nozzle, 

w = w = 0 , V = v„„, dkidy = d ddy = 0 

The near wall regions have two major characteristics that distinguish them from the bulk flow’ 
[31,35]; 

i. steep non-linear gradient in the flow properties, and 

ii. levels of Reynolds number are sufficiently low for local isotropy of small scale 
turbulence to prevail. 

Hence, the basic turbulence model becomes inapplicable in these regions and alternate 
representation must be sought. In the present study wall function procedure [31,35] has been 
used. The method is based on one-dimensional fluid flow, and assumes that the velocity profile 
in the near wall region is given by the universal logarithmic wall law which can be represented 
as : 


— = lln(y'R) 

M K 


(2.23) 


where u is the velocity parallel to the wall and { rj p)) is the friction velocity, x is the 

van Karmen constant and E is the roughness parameter typically taken as 9 for hydraulicalh- 
smooth walls, y~ is a dimensionless wall distance and is defined as : 


y 

11 


(2-24) 


where v is the normal distance from the wall. This law has been applied to a point whose v 
value is greater than 1 1.63. This is schematically shown in Fig.(2.2) . In the /^region specified 
above the Reynolds stresses are constant and the wall stresses, can be calculated as : 

(2.25) 




Wall 


Fig. (2.3) : Schematic representation of procedure aplplied for the near wall 
region on flow variable and turbulence parameter. 
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Further, in the near wall region the convection-diffiision of {um) is negligible. This implies that 
the rate of production of turbulent kinetic energy^ is equal to its dissipation. Such considerations 
lead to: 


_k_ 1 


(2.26) 


This relation is normally used as an boundary condition for k in the two equation k-e model. 
Similarly the boundary condition for € can be expressed as ; 


e = 



KV" 


(2.27) 


The two preceding relationships are to be applied to a near wall node whose y* is greater than 
11.63. 


For concentration field calculation zero flux boundary condition is imposed on all the solid 
surfaces, symmetrical planes and the free surface. In addition to this, .the following condition, 
i.e., 


at r = 0 , C,- = C° 


was applied to Eq.(2.22) as an initial condition, where, C“ is tracer concentration (tracer 
having the same density as the bulk liquid was considered) in the system at t=0. It is instructive 
to note here that one way coupling between Eq.(2.22) and the flow and turbulence model 
equations allows flow computation to precede the calculation of tracer dispersion.Using the 
converged flow field, the C, ( x,y,z,t) field was predicted. From the predicted concentiation at 
the outlet as a function of time, t, the C-curve was constructed. 

2.5 The General Differential Equation 

The general structure of the relevant differential equations ( Eqs.(2.1) to (2.6)) describing the 
conservation of mass and momentum and that of the k~ e Model of turbulence indicate that all 
the dependent variables of interest obey a generalised conservation principle. If the general 
variable is denoted by (p, the general differential equation can be written as; 



( 2 . 28 ) 


— (p4)) * div (pw<j)) = div (tgrad^) ^ 

The four terms in the general differential equation, are the transient term, convection term, 
the diffusion term (F is the diffusion coefficient ) and the source term, S^. The quantities 
/"and S are specific to a particular meaning of (p. . Depending on the general variable (p, 
an appropriate meaning will have to be given to the difhision coefficient Tand the source 
term S. For the present set of governing equations this is summarised in Table-2.3. 


Table-2.3 : Definition of T, (j) and S of Eq.{2.28). 


Conservation of 


F 


Mass 

1 

0 

0 

Momentum 

w 

Pr 

Eq.(2.2.a) 


mm 

Pe 

Eq.(2.3.a) 



P. 

Eq.(2.4.a) 

Turbulent kinetic energy 

k 

Pc/Ot 

Eq.(2.5.a) 

Dissipation rate of kinetic energy 

6 

Pe/0£ 

Eq.(2.6.a) 

Concentration 

m 

mmm 

Eq.(2.22) 


Hence, the procedure of casting any particular differential equation into the general form 
Eq.(2.28) requires manipulation of the equation so that the convection and diffusion terms 
conform to the standard form. The coefficient of grad cp is taken as the expression for F. 
and the remaining terms on the right hand side are collectively defined as the source term. 
S. Since, all the equations for mass, momentum transport: turbulence related phenomena 
and species conservation can be thought of a particular case of general equation Eq.(2.28) 
consequently, we need to concern ourselves with the numerical solution of Eq.(2.28) only. 

2.6 The Grid Arrangement 

A portion of the three dimensional grid is shown in Fig. (2.4). For the grid point P, points 
FandfVi east and west respectively) are its x-direction neighbours, while 2/ and S (north 
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Fig. (2.5) ; Scala.' and staggered control volumes in three dimensional Cartesian 
co-ordinate system 
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and south respectively) are the y-direction neighbours and T and B (top and bottom 
respectively) are the neighbours in the z-direction. The whole flow domain has been 
divided into a number of such non-overlapping control volumes. Staggered control volumes 
have been considered for the velocity components u, v and w. Fig.(2.5) represents the 
different staggered control volumes at which the discrete values of the velocities are 
located. These values represent the averages over the respective control volumes. Scalar 
variables i.e., p, k,ju^,p and e are stored at the grid nodes, while the velocities u, v and h’ 
are stored midway between the main grid nodes. The velocities are located where they are 
needed for the calculation of the convective fluxes of scalar variables. 

2.7 The Numerical Solution Procedure 

A control volume based finite difference procedure [34] has been adopted in the present 
study to solve the governing equations numerically. The solution is initiated by descritizing 
the domain in a large number of non-overlapping control volumes. The governing equation 
is integrated around all such control volumes to yield a system of algebraic equations, 
called the discretization equation [34], The governing equations Eqs.(2.1) to (2.8) have 
both first and second order derivatives. These terms are numerically integrated using the 
concept as applied to a combined convection-diffusion problem incorporating the hybrid 
differencing scheme proposed by Patankar [34]. The pressure-velocity interlinkage in 
Eqs.(2.2) to (2.4) have been solved by an implicit finite difference procedure referred to as 
SIMPLE [34] (Semi-Implicit Method for pressure linked equations). A typical discretization 
equation for the general differential equation (2.29) can be written as ; 


Ap 4)p = + Ajf <j)ppr + .4^^ + A^ 4)5 *Aj,4>t * Ag 4)^ + Ap 4>p + .vol (2.29) 

in which, Ap , is the centre point coefficient while Ap,A^, dy , As, Ap and Ap are the 
neighbour point ( three dimensional ) coefficients, representing the combined effect of 
fluid convection and diffusion. Ap°(t>p is zero for steady state problems. Further, Sq . 
represents the constant component of the linearized source term (viz., S = S^. + Sp (pp ] 
while the slope has been embodied in the expression of the centre point coefficient 



The centre point coefficient is represented as : 


4-4 ^ 4 . 4 . 4. 4 .A^.a;-S^.voI (2.30) 

Thus, in a system of n internal control volumes, n number of similar algebraic equations 
( viz., Eq.(2.30) are obtained. Further, since the discretization equations are obtained from 
the same governing equation Eq.(2.28), the former therefore embodies the same 
conservation principle as the latter one. This is an important feature of the control volume 
based formulation as against the routine Taylor series numerical solution procedure. 

Prior to the solution of the discretization equations the boundary conditions are also 
transformed into equivalent numerical form. Boundary conditions are only applied to 
control volumes located at the domain boundaries. Discretization equations are derived by 
the same procedure as mentioned. In essence, the incorporation of a boundary condition 
involves the modification of Ap and/or the source terms of the discretization equations of 
the boundary control volumes . 

The resultant set of discretization equations for both boundary and internal nodes were 
solved using the Tri-Diagonal Matrix Algorithm (TDMA) adopting a line by line solution 
procedure. In this, a particular grid line, say in z-direction, is chosen and assuming the 
dependent variable (^) to be guessed in thex and^ directions, the problem is thus reduced 
to a pseudo-one-dimensional situation and subsequently solved by TDMA. This was applied 
to all the grid lines in one direction and the entire process was repeated for the other tu’o 
space directions to obtain a tentative solution of the dependent variable , (p . This typical!}' 
constituted one iteration. The total number of iterations required was decided by the 
convergence criteria adopted, which in the present study was defined according to : 

E E E [Mp - ( ^ )] < (2.31) 


in which. 
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( 2 . 32 ) 


The triple sum in Eq.{2.31) represents the summation over the entire three dimensional 
domain. Under relaxation of dependent variable was employed to achieve/enhance 
convergence.If (p^ is the value from the previous iteration of the dependent variable (pp ( pp 
implies p evaluated at a point P) and a is the under relaxation factor then the change in the 
value of pp can be slowed down from iteration to iteration by the following mathematical 
treatment [34] ; 


j , cb , + 6 

= <[>;.« V ) (2.33) 

p 

in which, the value of p in the previous iteration. 

Under relaxation factors for the problem chosen were determined by trial and error. 
Typically, under relaxation factor of 0.5 was applied to ii, v and w and while a value of 0,7 
was used for the turbulence parameters, ^and e. No under relaxation factor was used for 
the calculation of concentration field. 

2.8 Modelling of the inclined tundish wall and flow modification devices 

To this end, the well laiown cell porosity technique of Moult and coworkers[36] have been 
applied. In this, the numerical calculation domain is considered to be perfectly rectangular 
and therefore, the grid is laid assuming the tundish to be rectangular Fig.2.6 (a). 
Subsequently, the inclined walls ( e.g., the exact tapered tundish geometry) are 
superimposed on the numerical calculation domain, so as to producer the true flow domain 
Fig2.6,(b). These are shown schematically in Fig. (2. 6). There, as seen, the intersection of 
the inclined tundish walls with the numerical grid results into two different kinds of control 
volumes in the numerical calculation domain. Active control volumes and inactive control 



volumes. Inactive control volumes are completely out side the flow domain ( shown with 
a dot at the centre in Fig.2.6(b)), Active control volumes while are either partly or 
completely within the flow domain. For the latter type of control volumes (e.g., completely 
within the flow domain and thus, entirely filled with liquid), the full form of the 
discretisation equation ( which is the algebraic equivalent of the partial differential 
equations, written below, for example, for a general variable cj) ) is applicable, viz., 

However, for the control volumes which are partly within the flow domain and partly 
outside, the coefficients of the discretisation equation are modified in such a way that only 
the correct, available control surfaces are applied to the calculation of convective and 
diffusive fluxes . Referring to Fig.(2.7) , thus, a set of blockage ratios or cell porosiw 
factors can be estimated (from knowing the tundish inclination and the characteristics of 
the grid system ) and these can be subsequently applied to modify the coefficients of the 
discretisation equations Eq.(2.34) appropriately. For a typical three dimensional control 
volume partly in the flow domain and partly outside ( see Fig.2.7), one can therefore, define 
a set of six different blockage ratios for the six faces of the control volume. On the basis 
of such, the following set of modified coefficients, applicable to such control volumes, 
can be conveniently estimated, viz.., 

^E, new ^E, old [ ^ ” ^E ] 

■^W, new old [ ^ ] 

^N, new ^N, old t ^ 1 

■^S, new ^S, old [ ^ 3 

new old [^ “ ^B 1 

•^T, new -^T, old [ ^ ] 

In the set of above expressions, b's are the blockage ratios in different directions ( as 
characterised by the subscript). 

An interesting situation arises, when all the blockage ratios for a given control volume 
assume a value of unity. This implies that the control volume is essentially outside the 



liquid liquid 

Fig. (2.6) : The characteristics of calculation and flow domains in terms of 
the finite difference grid. 

(a) calculation domain (b) flow domain 
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Fig. (2.7) : A control volume cut by the vessel wall (2-D representation). 


36 



flow domain and therefore inactive, as far as its influence on the predicted results is 
concerned . Under such a condition, the control volume or the characteristic grid node is 
completely isolated from the calculation domain by re-defming the components of the 
source terms, Sp and S, as follows: 

Sc =10^“ * the desired value 
and, Sp = -10^° 

Through such a practice, zero value of the dependent variables can be fixed at all the 
inactive control volumes. Similar approaches were adopted to model the flow control 
devices ( dams, weirs etc.) which are essentially inactive portion of the flow domain. 
Standard wall function procedures were applied on both the faces of the flow control 
devices which are normal to the longitudinal direction, x. 

2.9 The computer program 

The concept of casting all governing differential equations to a general differential 
equation has been outlined Section-2.5 . Consequently, a general purpose computer 
program was developed to solve Eq.(2.28). The flow diagram of the three dimensional 
model developed is shown in Fig.(2.8). 

The basic structure of the present combined Langrangian-Eulerian procedure has been: 

i. The total flow is equally divided between a preconsidered number of trajectories. 
On the basis of such, flow rate ( G) per trajectory is estimated. 

ii. The equations for bubble motion are solved for each of the trajectories and then 
superimposed on the Eulerian grid. From these, overall distributions of ai ,F^, 
and F^ etc. estimated. 

iii. Estimated , F, , F, and F, are incorporated into the equations of the continuous 
phase ( viz., the liquid ) and solved iteratively by embodying the SIMPLE procedure 
ofPatankar and Spalding [34]. 



iv. At the end of a predetermined number of iterations of the liquid phase equations, the 
trajectories of the bubbles are recalculated using the prevalent flow and turbulence 
fields. 

V. Based on this the distribution of a, , Fy and F^ etc, were updated. The above 
steps are repeated till the converged velocity field and bubble trajectories are 
obtained. 

vi. Using converged velocity fields, the concentration is calculated as a function of 
time. 

For single phase calculations ( viz. no gas injection), steps (iii) through (vi) were onl>' 
cosidered assigning a value of unity to a, . 



Fig. (2,8) : The flow diagram of the computer model 
























CHAPTER 3 


RESULTS AND DISCUSSION 

In the previous chapter, the salient features of the two phase, three dimensional turbulent 
flow model were presented. This chapter is concerned with the validation of the three 
dimensional two phase model and application of the model to fluid flow and Residence 
Time Distributions (RTD) in different tundish designs. 

3.1 Assessment of the computational prcedure developed 

The adequacy and appropriateness of the steady state, three dimensional, two phase 
turbulent flow model was assessed by carrying out some standard three dimensional flow- 
simulation tests. This included, (i) flows in enclosed cavity with one of the walls moving, 
(ii) entrance length v. Reynolds number calculation in ducts of square cross-section and (iii) 
prediction of bubble rise in a stagnant fluid. These calculations were t carried out under 
laminar flow condition ( by setting m- Me and modifying S„ , S,. and S,, appropriately in 
the flow equations). These are summarised below in detail. 

3.1.1 Laminar flow in an enclosed cubic cavity 

Steady, laminar , single phase flow in three dimension can be conveniently represented 
by the Navier-Stokes equations presented already. Thus, considering p, = \i, a i =1.0 and 
S^i = Sv = = 0 in the calculation procedure, fluid flow in a cubic cavity (I = Im, p = 

1 kg/m^ and p = 0.01 kg/m.s ) was predicted numerically considering one of the two 
horizontal cavity walls moving with a velocity of Im/s ( e.g.. Re=100) in the positive x 
direction . The boundary conditions for the test problem included, a zero velocity (no slip) 
on all the stationary walls. However, the influence of the mo-ving wall was incorporated 
via a viscous drag acting on the fluid in contact with the moving wall. In this one of the 
walls of the cubic cavity was made to move with a velocit}’ of \mj's in the positive x 
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direction as shown in Fig.3.1. The remaining five solid walls were considered to be 
stationary. In order to define the boundary conditions for the given problem exactly, a zero 
velocity (no slip) condition was applied on the stationary walls. In addition to this, to 
incorporate the influence of the moving north wall on the fluid flow, the general 
descritization equation for the dependent variable ( the x component of motion), i.e.. 

"p ’ “r "v "s "r "s ^ (3.1) 


was modified for the boundary control volumes lying in the immediate vicinity of the 
moving wall according to : 


■^p^p’° '^T^T (j.2) 

where, S\. = A„u„ + S,. and ii„ is the given velocity of the moving wall( = 1 m/s). 

The predicted flow on the central vertical XY plane (at z=0,5 m), transverse vertical YZ 
plane(at x=0.5 m) and transverse horizontal plane at (y=0.5) have been shown in Figs 3.2a 
through 3.2c. These flow fieldsare in good qualitative agreement with the equivalent 
computational results of Hwar and coworkers[37] . 

To illustrate this further, the predicted variations of the horizontal velocity component, u, 
along the central vertical line and that of the vertical velocity component along the 
horizontal centre line are shovra in Fig 3.3 There, equivalent results reported by Hwar and 
coworkers [37] have also superimposed for the sake of a direct comparison. Near perfect 
agreement between present and earlier estimates is readily apparent in Fig.3.3. It is 
important to mention here that results presented in Figs. 3.2 and 3.3. were derived using a 
31x31x31 uniform grid system which was essentially found to produce a grid independent 


solution. 



North 



X 


Fig. (3.1) : Schematic of the cubic cavity with one of the walls moving. 
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I.O m s 


Fig. (3.2a) : Predicted velocity field in longitudinal vertical plane at zAV of 0.5 
for cubic cavity test (Re=100) 
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(b) 


Fig. (3.2b) : Predicted velocity field for cubic cavity t3St (Re=100) 

(a) transverse vertical plane at x/L=0.5 

(b) tmsverse horizontal plane at y/H =0.5 





Dimensionless vertical velocity Dimensionless vertical distance 



(a) 



(b) 

Fig. (3.3) : Cubic cavity velocity profile for (Re=100) 

(a) vertical centre line 

(b) horizontal centre line 





3.1.2 Entrance length prediction for laminar flow 


When a flow enters a duct ( internal flow bounded by walls ) the viscous effects grow and 
boundary layer on the duct wall develops. The boundary layers eventually meet along the 
centreline and from this point onwards, the thickness of the boundary layers remain 
invariant in the streamwise direction. This is shown in Fig.(3.4) [38]. As seen, there is an 
entrance region where a nearly inviscid upstream flow converges and enters the tube. 
Viscous boundaiy' layers grow downstream, retarding the axial flow ii{x,r) at the wall and 
thereby, accelerating the centre core flow to maintain the incompressible, continuity 
requirement. At a finite distance from the entrance, the boundary layers merge and the 
inviscid core disappears. The tube flow then is entirely viscous, and the axial velociP’ 
adjusts slightly further until at .t = it no longer changes with x (no acceleration) and is 
said to be fully developed, i.e. u = u{r) only. Downstream of :c = the velocity profile as 
well as the wall shear stress is constant, and the pressure drop varies linearly with x, for 
both laminar and turbulent flows ( see Fig.3.4). 

Dimensional analysis [38] shows that the Reynolds number is the only parameter affecting 
entrance length. Thus, 

(3-3) 

D 

For laminar flow [39], the accepted correlation for entrance length in ducts is : 

^ = KRe (3.4) 

D 

Where A: is a constant of proportionality, the value of which, depends on the cross-section 
of the duct, e.g. for ducts with circular cross section, lc = 0.06. In the present work entrance 
lengths were calculated for ducts with square cross-section Fig.(3.5) for laminar flov 
condition. 

The dimension of square cross-section (D) was varied between 1x10 to 3x10 m 



Fig .(3 .4) : 



Schematic of the developing velocity profiles and pressure changes in the 
entrance of a duct flow. 
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(dimensions in the y and the z direction ). Dimension of the duct in the x direction (the 
direction of main flow) was defined such that x » 4 . Zero velocity (no slip) boundary 
conditions were specified for north, south, top and the bottom walls. Inflow and outfiou' 
boundary treatments were specified at the inlet (west) and outlet (east) planes respectively. 
Fig. 3.6 represents a. two dimensional grid arrangement near the entry and the exit planes. 
Since, u velocity (x component of motion) is known at the inlet plane, consequently, a 
numerical treatment analogous to that presented in Section-3 .1 A can be applied to 
accommodate the exact inflow boundary condition into the calculation procedure. For all 
the grid points P lying adjacent to the outflow boundary, the coefficient Ae is set to zero 
i.e. Peclet number is made sufficiently large Ae -0) to make the behaviour localh' 

one way in this region. In addition to this, the average exit velocity is specified at the 
outflow boundary on the basis of overall continuity which is updated from iteration to 
iteration by setting ii(mj,k) = ii{ni-l,j,k). 

Reynolds numbers (e.g., inlet velocity ) were varied arbitrarily between 25 to 180 and 
the flow field within the duct was predicted as a function of inlet Reynolds number. From 
the predicted results, entrance lengths were estimated graphically by plotting the centreline 
velocity as a function of x. Entrance lengths thus estimated were non-dimensionalised (by 
dividing L, by D) and plotted against the inlet Reynolds number. This is shown in 
Fig.(3.7). As shown a straight line was obtained between the two, satisfying the 
requirement dictated by in Eq.(3.4) (e.g., LJ Re). 

3.1.3 Prediction of bwbble rise in a stagnant fluid 

Assigning zero fluid velocity and turbulence in flow and turbulence model equations and 
considering a classical drag coefficient value of 8/3, the rise velocities of four different 
sized spherical cap air bubbles in water were estimated numerically following the 
computational procedure outlined already. An arbitrary vertical velocity was applied to 
Eq.(14) as the initial condition. The results as one would normally anticipate, indicated that 
the bubbles attain their terminal rise velocity within a short distance from above the point 
of injection. A comparison between numerically predicted terminal rise velocity and those 




© Scalar nodes 
— ► ’u’ velocity 

component 
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X 



Outflow boundary treatment 

Fig. (3.6) : The procedure adopted for modelling the inflow and outflow boundary 
conditions for the prediction of entrance lengths. 
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obtained on the basis of the following expression[32]. 


?7j,= 0.711. 


N P; 


(3.5) 


is illustrated in Fig.3.8. There, perfect agreement between the analytical and numerical 
estimates is at once apparent. Results presented so far evidently demonstrate the adequac^■ 
and appropriateness of the calculation procedure outlined in the preceding section. 


In addition to the aformentioned three tests, some other numerical experimentation s were 
also earned out ( viz., turbulent flow over a flat plate, cubic cavity with two parallel moving 
walls [40], These have also appear to confirm the reliability of the present computational 
procedure developed. 


3.2 Modelling of flow and Residence Time Distributions in steelmaking tundisb 
system 


The two phase, three dimensional turbulent flow model was applied to investigate flow and 
RTD in three different tundish designs. These included, (i) a single strand slab casting 
tundish, (ii) a two strand slab casting tundish and (iii) a six strand delta shaped billet 
casting tundish respectively. These are summarised below in detail. 


3.2.1 Flow and RTD in a single strand tundish 

As a starting case, fluid flew and RTD distribution behaviour in a water model tundish 
with single inlet and outlet and no flow control devices has been analysed computationalh'. 
The side walls of the tundish were tapered while the frontal wulls were perfectly vertical. 
The characteristic design features of the model tundish are summarised in Table 3.1 .A 
schematic of the tundish is shovra in Fig. (3.9). 
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Table-3.1: Characteristic parameters of the single strand model tundish [8] 


Parameter 

Numerical value 

Fluid 

Water, (//=0.001 kg m''s'‘,p=l000 
kglm^) 

Base length, m 

1.0 

Width (at base), m 

0.30 

Height, m 

0.37 

Liquid depth, m 

0.26 

Shape 

rectangular with sloping walls ( 1 5“) 

No. of strands 

one 

Ladle shroud diameter, m 

0.022 

Volumetric flow rate, mUs 

1.55 X 10^ 

Inlet Reynolds number 

9.0x10^ 


Computations were performed using a nonuniform grid covering the whole tundish. A grid 
system of 38 x 13 x 13 was chosen for the x (longitudinal), y (vertical) and the j 
( transverse) directions respectively. Under-relaxation factors of 0.5 for velocities and 0." 
for turbulence quantities were adopted . Details of the under-relaxation procedure employed 
has already been summarised in Chapter 2. A typical execution took almost 10,000 
iteration and this required about 300 min of CPU for flow field to converge 
Approximately, 150 min of CPU time was involved for the concentration to reach steady 
■ state level. All computations were carried out in the HP-9QQQ series machine at ITT. 
Kanpur. 

The prdicted three dimensional flow field in the single strand tundish is shown in Figs. 3.10 
through (3.12) via a series of two dimensional plots. Figure (3.10) shows fl.ow at different 





Fig.(3.9) : Schematic of the single strand model tundish dimensions[8]. 


elevations in the transverse horizontal planes (X-Z planes). There the bottom most figure 
(3.10d) corresponds to the flow on a horizontal plane that is closese to the tundish base, 
while the upper most plot (3 . 1 Oa) corresponds to the flow in the vicinity of the free surface. 
Referring to Fig.3. 10, it is seen that the jet following its impingment on to the tundish base 
spreads radially along all directions. Further more, while the flow closer to the vertical 
frontal wall is directed primarily towards the exit nozzle, in contrast shows a flow reversal 
along the central line of the particular plane. In the vicinity of the exit nozzle, a sink like 
vortex can be observed which is due to drainage effect. As one moves towards the free 
surface, this flow characteristic ( flow pattern on plane (a)) changes and more liquid is seen 
approaching the input stream from towards the right of the vessel. Finally, near the free 
surface flow pattern again shifts where surface flow primarily directed towards the exit 
nozzle is observed. It is important to note here that predicted flow on the horizontal planes 
exhibit symmetry with respect about the horizontal centreline. 

In Fig. 3.11 predicted flows on four different vertical planes are shown, because the tundish 
is symmetrical about the central verical plane consequently flow on the frontal half has 
been shown in Fig.3. 1 1 . On the central vertical plane Fig.3. 1 la, it is seen that the jet upon 
entry, entrains surrounding fluid and undergoes considerable expansion before impinging 
on to the base of the tundish. Towards the left of the incoming jet, a disticnt recirculatoiy 
loop is seen ( created by the confmemet of the jet by the tundish vertical and side walls). 
Towards the bottom of the same plane, the flow is seen to be first directed towards the 
orifice , However at some horizontal distance from the inlet plane, the flow is seen to be 
reversed. The flow pattern changes drastically as one visits various planes from the 
centreline to the proximity of the frontal vertical wall. In Figs. 3. 11(b) and 3.11(c) 
considerable amount of fluid is seen to be moving towards the jet impinging area. Further 
more, the surface directed flow becomes more intense near the frontal vertical wall ( e.g., 
as one moves down from Fig.3.1 1(a) to 3,1 1(d). Similarly, in Fig. 3.12, flow patterns along 
the transverse vertical planes are shown. This, shows the distinct recirculatoiy loop created 
by the confinement of the jet between the two frontal walls particularly, in the proximift' 
of the input stream. 


















0.4 m/s 


- — -lx +iTnHi<5Vi in the transvese plane at 

Fig.{3.12) : Predicted flow field in (b) 0.28 (c) 0.67 (d) 0.78 

dimensionless distance (x/L) oi ( J — 
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As pointed out already, the Residence Time Distribution characteristics in the mode! 
tundish was investigated through numerical solution of the convection- turbulent diffusion 
equation summarised already in Chapter 2. In this, at time t=0 , the inlet control volume 
was made full of tracer of the same density as that of bulk liquid phase (e.g., water) and 
the concentration of the tracer was then computationally monitored as a function of time 
at the control volume passing through the exit noz 2 de and one control volume vertically 
above it. The concentration of tracer at these these locations was subsequently non- 
dimensionalised using the average concentration that the tracer would assume in the vessel, 
if the later were perfectly mixed in the absence of any drainout ( = {mass, Volume 
) = 0.02 kg/m^). Similarly, time was also non-dimensionalised using the nominal holding 
time Eq.(l . 1) of the fluid in the vessel ( /;). The result thus obtained ( C/C„ vs. t/ Q was 
then plotted to get the desired C-curves. In Fig.(3.13), the C-curves obtained 
computationally for the control volume just above the exit nozzle is shown. There as seen 
the tracer is initially of negligible concentration at the exit stream, but it quickly increases 
to a maximum, before showing a gradual fade out. These features of the RTD curve is a 
characteristic of the metallurgical tundish systems and is similar to those reported by 
numerous other investigators [8,12,18]. 

The minimum breakthrough time, the time at which peak concentration is reached and the 
mean residence time were calculated from the "C" curve. Mean residence time has been 
calculated using eq(1.2). On the other hand, minimum break through time and Peak 
concentration time were deduced from visual examination of the predicted C curve. 
Residence Time Distribution characteristics thus obtained have been compared directly in 
Table 3.2 which also contains corresj 3 onding experimental measurements reported by Singh 
and Koria. There, excellent agreement between prediction and experiment is readily 


apparent. 



Dimensionless concentration 





Table 3.2 : Chracteristics of the RTD in the single strand tundish and their 

comparison with equivalent measurements reported in the literature[8]. 



3.2 Flow and RTD in a two strand slab casting tundish 


Table 5.5: Characteristic parameters of the two strand model tundish (13) 


Parameter 

Numerical values 

Fluid 

Water, (//=0.001 kg m'^s'^ ,p=\000 
kglm^) 

Base length, m 

1.2 

Width (at top), m 

0.2 

Width at the base ,m 

0.154 

Liquid depth, m 

0.14 

Shape 

rectangular ( -10° inclined) 

No. of strands 

two 

Tundish nozzle diameter, mm 

10 

Volumetric flow rate, m^/s 

0.93 X 10-^ 

Inlet Reynold's Number 

9.43 X 10^ 












Three different configurations of the two strand tundish were used for numerical 
computations. In each configuration, the operating parameters were kept identical to those 
summarised in Table 3.3. Schematics of the three configuration are shown in Fig.3.14. 
Owing to symmetry, computations were restricted to only one quarter of the tundish. A 
32x12x13 grid system was applied to seek the numerical solution. Details of the computaed 
resuls obtained together with discussion are summarised below. 

In Figs.3.15 through 3.16 flow fields in the absence of any flow modification devices are 
shown through a series of two dimensional plots. In Fig.3. 15, flow fields are presented for 
three different vertical planes. There, the the central vertical plane shows strong downward 
fluid motion in the vicinity of the jet entry location. In the remaining of the plane, fluid 
is relatively stagnant gaining momentum again in the proximity of the exit nozzle. 
Furthennore, the flow velocities particularly towards the bottom of the tundish are seen to 
be alligned in a direction towards tire exit nozzle for some distance only. However, as one 
approaches the frontal wall, the intensity of the exit nozzle directed flow increases 
substantially Fig.3, 15(c), These chracteristics of the flow are qualitatively very similar to 
those presented earlier for the single strand tundish . This is to be expected since there is 
no significant geometrical difference between a single strand and one half of a two strand 
tundish in a qualitative sense. 

In Figs. 3.16(a) through (c), flow pattern in three different horizontal planes have been 
shown for the operating conditions summarised in Table 3.3.. This shows that the jet 
following its im oingement on the tundish base, flows outwards . A portion of the jet flows 
towards the frontal wall and then rises up ( see also Fig.3. 15(c)) while the remaining 
portion of the j(h flows preferentially towards the exit nozzle. As one moves towards the 
free surface this chracteristics of the flow pattern changes substantially and ultimately near 
the free surface flow is seen to be directed towards the jet entry region. In the horizontal 
planes, except for the immediate vicinity of the inlet stream, flow is seen to be relatively 
sluggish. 






Fig.(3.14) : 


o U .f three different configurations of the tNvo strand tundish used for 

simulation (a) without baffle (b) with a dam (c) with a dam and a wtet 

(c) with a bubbler in place of dani 



In Fig.3. 1 8, flow field in three different vertical planes are shown for the two strand tundish 
with a dam placed at a distance of 0.36m from the central plane of symmetry . These figures 
show that the flow pattern, particularly between the inlet stream and the dam is not much 
affected by the presence of the dam within the flow domain. However in the vicinity of the 
dam and beyond, the flow pattern appears to have changed with respect to an equivalent 
unbaffled situations. Similarly, in Fig.3. 1 8(a) through (c) flow patterns in three different 
horizontal planes are shown. There, flow fields on those planes above the dam are 
practicaly identical to those presented earlier for an equivalent no-baffle situation. The dam 
seems to have no influence on the flow patterns generated on the horizontal planes 
particularly in regions bounded by the the inlet steam and the dam. This is to be expected 
since the dam is located much down stream of the inlet nozzlle. 

In Figs. (3.19) through (3.20) flow fields in the two strand tundish with a dam and a weir 
is shown. The weir as seen from these figures appear to influence the predicted results 
considerably because of its proximity to the jet entry location. Interestingly, large flow 
velocities coupled with appreciable fluid turbulence is seen to be primarily restricted 
between the region bounded by the inlet stream and the weir, which is desirable feature as 
far as tundish hydrodynamics are concerned. 

Finally, flow fields were computed in the two strand tundish with gas ( Q= 2 lit/ min) 
stirring through a plug located at the base of the vessel. These were carried out in the 
absence of any flow modification devices. Predicted flow patterns are shown in Figs 3 .21 
through 3.22. There the general features of the flow is very similar to those presented 
earlier for an equivalent no baffle situation. Only in the proximity of the bubbler, some 
noticable changes in the flow intensity can be seen and bubbler forms recirculatory motions 
on its either sides. 

Corresponding predicted Residence Time Distribution chracteristics ( see Fig.3. 23) are 
shown in Table 3.4 for the four different cases . There, predicted minimum break through 
time, time at which peak concentration is reached and the average residence time have been 
summarised. It is clear that the dam and the dam+weir configuration tends to increase the 






Fig.(3.15) : 


Prdicted flow field in the two 
longitudinal vertical planes at 


strand symmetrical tundish without baffle in the 
dimensionless distance (zAV) of (a) 0.99 (b) 0.70 


(c) 0.34 
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Fig(3.16) 


(c) 


0.011 m/s 


Predicted flow field m the too strand s^etrictd ^ , 

horizontal transverse planes at different dimensmless height (y/H) of (a) O.PD (b) 

0.50 (c) 0.05 
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(C) 


0.3 xnJs 


Fig.(3.17) : 


Prdicted flow field in the two strand symmetrical tundish with a ^ 
longitudinal vertical planes at dimensionless distance (z/W) of (a) 0.99 (bj U. 
(c)0.34 
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(C) 


0.011 m/s 


Fig(3.18) : 


(b) 0.49 (c) 0.05 






.99 (b) 0.70 (c) 0.34 




(a) 



(b) 
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0.011 m/s 


Fig. (3.20) ; Predicted flow field in the two strand symmetrical mndish with a wier and a dam in 
horizontal transverse planes at different dimensinless heights of (a) 0.95 (b) 0.49 
and (c) 0.05. 
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0.3 mJs 

(c) 


Fig.(3.20) : Prdicted flow field in the two strand symmetrical tundish with a bubbler in 

the longitudinal vertical planes at dimensionless distance (zAV) of (a) 0.99 (b ) 
0.70(c) 0.34 



(b) 



(c) 


0.011 m/s 


Fig(3.21) ; 


Predicted flow field in .he two st^d W »”« 

the horizontal transverse planes at different dimensinless heignt ty/ ) 

(b) 0.50 (c) 0.05 





Dimensionless concentration Dimensionless concentration 


With a dam and a wier 
With a dam 
Without baffle 




Fig.(3.22) 


(b) 

RTD for two strand tundish for four different case (a) wthoiri baffle (b) with a 
to (oT2a weir and a dam (d) with abubbler located m place of a dam 




residence time of the fluid in the tundish as one would normally anticipate. Significant!}-, 
the bubbler seem to serve no purpose as far as increase of residence time of fluid in the 
tundish are concerned. In fact gas stirring needs to be applied in conjunction with flow 
modification devices to obtain the necessary benefits. This is however left for a future 
study. 

For the three configurations the dead volume, dispersed plug volume and complete mixina 
volumes were calculated from the eqs. 1.3 to 1.5 and compared with the equivalent 
experimental values and a good agreement can be observed. 


Table 3.4 : Residence Time distribution characteristics in the two strand tundish 
as a function of various configurations 


Configuration of the 

tundish 

RTD characteristics 

Numerical values, 

s 

Without baffle 

Minimum breah through time 

5 

Peak concentration time 

no 

Mean residence time 

200 

with a dam only 

Minimum breah through time 

14 

Peak concentration time 

115 

Mean residence time 

213 

with a weir and a dam 

Minimum breah through time 

14 

Peak concentration time 

115 

Mean residence time 

211 

with a gas bubbler 

Minimum breah through time 

06 

Peak concentration time 

40 

Mean residence time 

141 


141 
































Table 3.5 : Comparision of flow characteristics of two strand tundish with 
experimentally observed values. 


Configuration of 

Flow characteristics 

Predicted 

Experimenta 

the tundish 


values (%) 

-1 values (%) 

Without baffle 

dead volume 

26 

29 


dispersed plug volume 

19 

3 


complete mixing volume 

55 

68 

With a dam only 

dead volume 

21 

7 


dispersed plug volume 

23 

7 


complete mixing volume 

56 

86 

With a weir and a 

dead volume 

22 

15 

dam 

dispersed plug volume 

23 

8 


complete mixing volume 

55 

77 


3.3 Flow and RTD in the six strand billet casting tundish 

Finally flow and RTD simulation were carried out in a delta shaped six strand billet casting 
water model tundish. Schematics of the six strand tundish is shown in Fig. 3.24. The 
operating chracteristics of the six strand tundish are summarised in Table 3.5. 
Computations were perlbrmed using a nonuniform grid covering the the half of the tundish 
because of symmetry . The grid set of 40 x 15 x 20 were chosen for the x (longitudinal), 
y (vertical) and the z (transverse) directions respectively. A tvpical execution involving 
8000 iteration required about 350 min of CPU time for flow field to converge and an 
additional 200 min of CRf/time for the concentration field to attain steady state. 



























Fig.(3.23) Schematic of the six strand tundish with all dimensions 


Tabie 3.6 : Characteristic parameter of the six strand delta shape water model 
tundish. 


Parameter 

Numerical values 

Fluid 

Water, (//=0.001 kgm:^s'^,p=\Q00 
kglm^) 

Base length, m 

2.86 

Width (at base), m 

0.145 at side walls and 0.225 at the 

centre 

Height, m 

0.330 

Liquid depth, m 

0.275 

Shape 

delta shape with sloping walls (10°) 

No. of strands 

six 

Ladle shroud diameter, m 

0.040 

exit nozzle dia,/« 

0.012 

Volumetric flow rate, n^/s 

6.0 xlO"* 


Detai .3 of the results and analyses on predicted flow and RTD are presented below. 

The velocity profiles near the frontal wall is predominantly tow^ds the outlet side which 
is due to large drainage effect associated with the three outlets ( Fig.3.25d). There some 
part of fluid is seen to be directed towards the free surface due to entrainment of the liquid 
by the down coming jet and this forms a recirculatoiy motion to the right. Fluid velocities 
on the plane passing through the outlets are very weak elsewhere except in the proximity 
of the outlets( Fig.3.25c) and a very weak recirculation cmi be observed close to bottom 




wall and at the inlet stream. As seen from Fig. {3.29b) entrainment has become quite 
strong as is evident from fluid flowing towards the inlet jet close to free surface in the 
right half of the tundish. Wliile in the left part of the tundish we can see fluid comes from 
the free surface to the bottom because of the drainage effect of the outlets. As seen in 
Fig.(3.25a), which has the inlet stream the fluid is stagnant everywhere except in the region 
of inlet Jet, where a strong recirculation can be observed because of combined entrainment 
effect of inlet jet and drainge effect of the outlets. 

In transverse horizontal planes, close to bottom Fig(3.26d) flow is predominantly towards 
the exit nozzles.The incoming jet strikes at the base of the tundish and spreads radially 
outwards. More fluid flows towards the frontal wall because of the drainage effect of the 
outlets. Fluid velocities are weak near exit nozzles # 2 and # 3 which are further from the 
flow inlet. As one moves from the bottom plane towards the free surface the reversal of the 
fluid flow reversal can be witnessed. In Fig.(3.26b) on the right of the tundish, in the 
middle of the tundish the flow is mainly towards the inlet jet because of the etrainment, 
whereas in the region close to the walls and to the left: of the tundish the flow is towards 
the exit side. As we move towards the free surface, to the left of the tundish the flow is 
predominantly towards the exit side and to the right the flow is mainly towards the inlet jet 
with high velcities, because of the large entrainment effect. 

Fig.(3.27), shows the RTD characteristics of the given configuration. Table (3.7) shows 
the comparision of minimum breakthrough time of prdicted and experimentaly observed 
values. Only minimum breakthrough for three outlets have been compared with the 
experimentally observed values, as only those experimental values were available by the 
time this study was completed. As we can see that we have excellent agreement of the two. 
As expected the minimum break through times for exit #1 is the shortest while the same for 
exit nozzle #3 is the largest. This is to be expected since the exit nozzle # 1 is closest while 
exit nozzle #3 is farthest from the inlet stream. 




(d) 


P^ediCd flow flew six 

vertical planes at dimensionless distance z/ w 01 w w 








(C) 

Fig. (3.25) Predicted flow field in six strand water model tundish ( no baffle) in transvene 

horizontal planes at dimensionless height y/H of (a) 0.05 ( 3 ) 0.50 (c) 0.95 


Dimensionless concentration 
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Comparison of RTD for six 
different outlets 


strand delta shaped tundish without baffle for 



Table 3.7 : Comparision of predicted and experimentally observed values of 

minimum breakthrough times for different outlets of six strand delta 
shaped water model tundish 


Exit nozzle 

1 

RTD characteristics 

i 

Prdiected values 

Experimental 

values 

1 

i 

Minimum break through 

time 

2 sec 

4 sec 

2 

Minimum break through 

time 

18 sec 

1 

20 esc 

3 

1 

Minimum break through 

time 

100 sec 

120 sec 














CHAPTER 4 


CONCLUDING REMARKS 

The following three major conclusions can be derived from the present study; 

(i) A steady state state, three dimensional, two phase, turbulent flow calculation 
procedure has been developed in Cartesian coordinate system adopting the control 
volume based finite difference procedure in conjunction with the SIMPLE algorithm of 
Patankar and Spalding. The calculation procedure developed has been tested against 
several bench mark solutions to essentially demonstrate its adequacy and appropriateness. 

(ii) The model has been modified through incorporation of cell porosity technique for 
flow and RTD simulation in three different designs of tundish which are most commonly 
encountered in the steelmaking industry. 

(iii) RTD predictions for all tundish geometries have been compared with equivalent 
experimental results published in the literature and good agreement has been achieved. 



CHAPTER 5 


RECOMMENDATIONS FOR FUTURE WORK 

The scope of the present work can be further expanded considering the following: 

(i) The two phase, three dimensional, turbulent flow model which has been developed in 
the present work can be suitably modified to carry out equivalent computations in polar 
co-ordinate system. The code can then be used to simulate a variety of metallurgical 
processes such as, asymmetric gas injection in vessels, side blown converters, R.H 
degassers and so on. 

(ii) The influence of temperature variations on fluid flow (free convection affects) has 
been neglected in the present study. Integration of heat transfer aspects to the present 
model will enable it to predict the affect of temperature losses or auxiliary heating on fluid 
flow. Similarly, appropriate form of Maxwell equations can be incorporated in the body 
force terms of the Navier-Stokes equation in the present model, to study the affect of 
magnetic field on fluid flow in metallurgical tundishes and its resultant influences on the 
tundish efficiency. 

(iii) The gas injection device is to be incorporated in conjunction with appropriate flow 
modification devices ( dam and weir). Instead of a single bubbler, a set of bubblers is to 
considered in the model to investigate the benefits of gas injection. 

(iv) Extensive numerical trials are to be carried out to study the influence of grid size, 
incremental time step and the distance of the near wall node on the predicted results. 
Similarly, under relaxation, solution of the discretisation equations etc. are to be looked 
at more rigorously to enhance the computational efficiency of the calculation procedure 



developed. 
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